mean_cost_diff_nalox <- mean(m_outcomes_psa[, "cost_diff_nalox"])
ci_cost_diff_nalox <- quantile(m_outcomes_psa[, "cost_diff_nalox"],
probs = c(0.025, 0.975, 0.5))
mean_cost_diff_ss <- mean(m_outcomes_psa[, "cost_diff_ss"])
ci_cost_diff_ss <- quantile(m_outcomes_psa[, "cost_diff_ss"],
probs = c(0.025, 0.975, 0.5))
mean_cost_diff_pg <- mean(m_outcomes_psa[, "cost_diff_pg"])
ci_cost_diff_pg <- quantile(m_outcomes_psa[, "cost_diff_pg"],
probs = c(0.025, 0.975, 0.5))
mean_cost_diff_all <- mean(m_outcomes_psa[, "cost_diff_all"])
ci_cost_diff_all <- quantile(m_outcomes_psa[, "cost_diff_all"],
probs = c(0.025, 0.975, 0.5))
mean_cost_diff <- c(paste0(round(ci_cost_diff_nalox[3]/1000000, 0), " [",
round(ci_cost_diff_nalox[1]/1000000, 0),", ",
round(ci_cost_diff_nalox[2]/1000000, 0), "]"),
paste0(round(ci_cost_diff_ss[3]/1000000, 0), " [",
round(ci_cost_diff_ss[1]/1000000, 0),", ",
round(ci_cost_diff_ss[2]/1000000, 0), "]"),
paste0(round(ci_cost_diff_pg[3]/1000000, 0), " [",
round(ci_cost_diff_pg[1]/1000000, 0),", ",
round(ci_cost_diff_pg[2]/1000000, 0), "]"),
paste0(round(ci_cost_diff_all[3]/1000000, 0), " [",
round(ci_cost_diff_all[1]/1000000, 0),", ",
round(ci_cost_diff_all[2]/1000000, 0), "]"))
mean_death_diff_nalox <- mean(m_outcomes_psa[, "death_diff_nalox"])
ci_death_diff_nalox <- quantile(m_outcomes_psa[, "death_diff_nalox"],
probs = c(0.025, 0.975, 0.5))
mean_death_diff_ss <- mean(m_outcomes_psa[, "death_diff_ss"])
ci_death_diff_ss <- quantile(m_outcomes_psa[, "death_diff_ss"],
probs = c(0.025, 0.975, 0.5))
mean_death_diff_pg <- mean(m_outcomes_psa[, "death_diff_pg"])
ci_death_diff_pg <- quantile(m_outcomes_psa[, "death_diff_pg"],
probs = c(0.025, 0.975, 0.5))
mean_death_diff_all <- mean(m_outcomes_psa[, "death_diff_all"])
ci_death_diff_all <- quantile(m_outcomes_psa[, "death_diff_all"],
probs = c(0.025, 0.975, 0.5))
mean_death_diff <- c(paste0(round(ci_death_diff_nalox[3], 0), " [",
round(ci_death_diff_nalox[1], 0),", ",
round(ci_death_diff_nalox[2], 0), "]"),
paste0(round(ci_death_diff_ss[3], 0), " [",
round(ci_death_diff_ss[1], 0),", ",
round(ci_death_diff_ss[2], 0), "]"),
paste0(round(ci_death_diff_pg[3], 0), " [",
round(ci_death_diff_pg[1], 0),", ",
round(ci_death_diff_pg[2], 0), "]"),
paste0(round(ci_death_diff_all[3], 0), " [",
round(ci_death_diff_all[1], 0),", ",
round(ci_death_diff_all[2], 0), "]"))
mean_oddeath_diff_nalox <- mean(m_outcomes_psa[, "oddeath_diff_nalox"])
ci_oddeath_diff_nalox <- quantile(m_outcomes_psa[, "oddeath_diff_nalox"],
probs = c(0.025, 0.975, 0.5))
mean_oddeath_diff_ss <- mean(m_outcomes_psa[, "oddeath_diff_ss"])
ci_oddeath_diff_ss <- quantile(m_outcomes_psa[, "oddeath_diff_ss"],
probs = c(0.025, 0.975, 0.5))
mean_oddeath_diff_pg <- mean(m_outcomes_psa[, "oddeath_diff_pg"])
ci_oddeath_diff_pg <- quantile(m_outcomes_psa[, "oddeath_diff_pg"],
probs = c(0.025, 0.975, 0.5))
mean_oddeath_diff_all <- mean(m_outcomes_psa[, "oddeath_diff_all"])
ci_oddeath_diff_all <- quantile(m_outcomes_psa[, "oddeath_diff_all"],
probs = c(0.025, 0.975, 0.5))
mean_oddeath_diff <- c(paste0(round(ci_oddeath_diff_nalox[3], 0), " [",
round(ci_oddeath_diff_nalox[1], 0),", ",
round(ci_oddeath_diff_nalox[2], 0), "]"),
paste0(round(ci_oddeath_diff_ss[3], 0), " [",
round(ci_oddeath_diff_ss[1], 0),", ",
round(ci_oddeath_diff_ss[2], 0), "]"),
paste0(round(ci_oddeath_diff_pg[3], 0), " [",
round(ci_oddeath_diff_pg[1], 0),", ",
round(ci_oddeath_diff_pg[2], 0), "]"),
paste0(round(ci_oddeath_diff_all[3], 0), " [",
round(ci_oddeath_diff_all[1], 0),", ",
round(ci_oddeath_diff_all[2], 0), "]"))
effects_tbl <- cbind.data.frame(mod_names[-1], mean_cost_diff,
mean_death_diff, mean_oddeath_diff)
effects_tbl |>
flextable() |>
add_header_row(colwidths = c(1, 3),
values = c("", "Change Compared to No Intervention \n Median [95% Credible Intervals]")) |>
set_header_labels(values = list(
`mod_names[-1]` = "Interventions",
mean_cost_diff = "Discounted Net Present \n Costs in Millions ",
mean_death_diff = "Deaths ",
mean_oddeath_diff = "Opioid-related Overdose Deaths")) |>
theme_vanilla() |>
set_table_properties(layout = "autofit") |>
align_text_col(align = "center", header = TRUE, footer = TRUE)